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Abstract. We present spectro-astrometric observations for 28 southern pre-main sequence (PMS) stars and 
investigate their circumsteUar environment down to AU scales. The structures detected in the "position spectra" 
include: (1) almost all the known binary companions in our sample (Sz 68, Sz 41, HO Lup, VW Cha, S CrA, AS 
205), (2) companion candidates which have not been detected by infrared speckle techniques (T CrA, MWC 300), 
(3) monopolar and bipolar jets (AS 353A, CS Cha), (4) a combination of jets and a bow shock (VV CrA), and 
(5) a combination of a jet and stellar companion (R CrA). Results in known binaries show that this technique is 
capable of detecting binaries with separations down to ~10 milliarcsec. Both components in each binary appear 
to have strikingly similar profiles in Ha emission, indicating a similarity of circumsteUar activity (mass accretion 
and /or a wind), and supporting the scenario of core fragmentation for the mechanism of binary formation. The 
bipolar Ha jet in CS Cha has a spatial scale of ~1.5 AU, similar to that previously observed in RU Lup, and likely 
be heated by a mechanism other than shocks. From the spatial scale, velocity, and Ha luminosity, we estimate 
the mean hydrogen density in the AU-scale bipolar flows to be <^10 7 cm -3 . The bipolar geometry in these jets 
can be explained by the presence of a disc gap/hole at AU scales, which could be induced by a gas-giant planet 
at the ice condensation radius. 

Key words, accretion, accretion discs - line: formation line: profiles stars: activity stars: evolution stars: pre- 
main-sequence 

1. Introduction out by photometry and spectroscopy without any direct 

spatial information, and optical-to-radio interferometry in 
The circumsteUar environment within 10 AU of young a unu t e d number of objects 
stars are of particular interest for star and planet forma- 
^ ■ tion. In this region, the star accretes the material from the To obtain s P atial information well below these resolu- 
circumstellar disc, and drives an energetic jet/wind (see, tl0ns ' we have be S un "spectro-astrometric" observations 
e.g., Najita et al. 2000; Shu et al. 2000, Konigl & Pudritz of P re " mam sequence (PMS) stars. Bailey (1998a) applied 
2000). In the circumsteUar disc, planets would be grow- this technique for a PMS star, evolved star, and AGN, 
ing up and tidally interact with the circumsteUar environ- showin S that the technique allows us to study the spa- 
ment (see e.g., Lubow & Artimowicz 2000). This region tial structure of their emission/absorption features down 
also includes a stellar companion in a significant num- to milliarcsec scales. Bailey (1998b) succeeded to detect 
ber of objects, and the population, mass distribution, and almost a11 known PMS bmanes m his sam P le ' P rovln S the 
orbital parameters of these companions could give use- Potential of this technique to find new binaries with sep- 

ful constraints on the understanding of binary formation arations of wel1 below °- 1 arcsea Takami et aL ( 2001 ) 

(see e.g., Mathieu et al. 2000). However, even the Hubble then revealed the P resence of bi P olar structure at -30 

Space Telescope or adaptive optics on 10-m telescopes mas scale m the outflow of the T Taun star RU Lup. This 

cannot resolve this region: a spatial scale of 10 AU corre- bl P olar structure contrasts with the jet at larger scales, 

sponds to 0" .07 in the nearest star forming regions, which which exhlblts onl y blueshifted components, suggestive of 

is comparable to the angular resolutions of these facilities. the Presence of a disc gap at a few AU which allows the 
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Observational studies in this region have thus been carried 



redshifted Ha outflow to be seen. 

In this paper, we present results of spectro-astrometric 

Send offprint requests to: Michihiro Takami, e-mail: observations for 28 southern PMS stars, and investigate 
takami@star.herts.ac.uk their binarity, outflow activity, and disc structure down 
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to AU scales. Section 2 summarizes the observations. 
Sections 3 shows the detected structures due to stellar 
companions, jets, and bow shocks. Section 4 shows the 
Ha profiles for each component of known binaries, which 
probes the similarity or difference of circumstellar activ- 
ities between binary components. Section 5 investigates 
the detectability of PMS binaries, the mechanism of bi- 
nary formation, physical conditions of the Ha outflows 
down to AU scales, and the possibility of planet forma- 
tion in a few PMS stars. The conclusions of this paper is 
described in Section 6. 



2. Observations 

Observations were carried out at the 3.9-m Anglo- 
Australian Telescope using the RGO spectrograph. The 
configuration with a 1-arcsec slit width, a 1200 line mm -1 
grating, and the 82-cm camera provided a spectral reso- 
lution (A/AA) of 7000. The pixel scale of 0.15 arcsec with 
the MITLL CCD (2048x4096 deep depletion CCD from 
MIT Lincoln Labs) provides good sampling of the see- 
ing profile, whose FWHM varies between 1 and 3 arcsec 
during the observations. The flat fields were made by com- 
bining many exposures with the spectrograph illuminated 
by a tungsten lamp. Wavelength calibration were made by 
observations of a CuAr lamp. 

Table 1 shows the log of the observations. The tar- 
gets were selected based on their brightness (V 5s 14), 
and most of them lie in the nearest low-mass star form- 
ing regions including Lupus (190 pc — Wichmann et 
al. 1998), Chamaeleon (170 pc — Bertout et al. 1999), 
Corona-Australis (130 pc — Marraco & Rydgren 1981), 
Ophiuchus (140 pc — de Geus & Burton 1991), and Upper 
Scorpius (150 pc — de Zeeuw et al. 1999). The spec- 
tra were obtained with a wavelength coverage of 6480- 
6780 A (but 6140-6800 A for VV CrA) to observe lines 
including Ha, [Sn], Hei 6678 A, and obtained at four 
slit position angles (0°, 90°, 180°, and 270°) to deter- 
mine the positional displacement in N-S and E-W direc- 
tions. Although several targets were previously observed 
by Bailey (1998b), we observed them again to ensure all 
the data had uniform quality with which to carry out a 
statistical study. 

The data were reduced using the FIGARO package 
provided by Starlink. After subtracting the bias level and 
flat-fielding, the position spectrum was determined by fit- 
ting the seeing profile at each wavelength with a Gaussian 
function. After the fitting, the position spectra with op- 
posite position angles (0°-180°, or 90°-270°) were sub- 
tracted to remove instrumental effects (see Bailey 1998a. b 
and Takami et al. 2001 for detail). Each position spectrum 
has an arbitrary zero point, which we adjust to correspond 
to the continuum position. In addition to the position 
spectra, the intensity spectra were obtained by subtract- 
ing the bias, flat-fielding, subtracting the adjacent sky and 
extracting bright columns on the CCD. The systemic ve- 
locity of the targets was calibrated using photospheric ab- 



sorption lines including Li I 6707.82 A, Nil 6613.62 A, Cai 
6439.08 A, Fei 6475.63 A, and Fen 6456.38 A. 

The measured positional accuracy, cr, depends on the 
seeing size and the photon number at each wavelength. 
Each photon should be within a distance of the order of 
the seeing size from the actual position of the source, and 
this positional uncertainty is reduced by a factor of N 1 / 2 
by observing N photons. Consequently, the detection limit 
is described as follows (Bailey 1998a): 

cr(milliarcsec) = const. ■ FWHM(milliarcsec) • iV~ 1//2 , (1) 

where FWHM is the full width half maximum of the see- 
ing. Table 1 provides the constant value of 0.5±0.2, in 
agreement with that estimated using a Gaussian seeing 
profile (0.4). These results imply that the positional accu- 
racies in our spectra are dominated by photon noise, other 
noise (e.g., readout noise) being negligible. 

3. Positional displacement 

Displacement in the position spectra of 12 objects were 
detected, including: (1) known binaries with separations of 
0.1-1.5 arcsec (Sz 68, Sz 41, HO Lup, VW Cha, S CrA, and 
AS 205); (2) T CrA and MWC 300, whose displacement 
can be explained by a stellar companion; (3) CS Cha, AS 
353A, and VV CrA, whose displacement can be attributed 
to outflowing gas; and (4) R CrA, whose displacement 
can be attributed to a combination of a stellar companion 
and outflowing gas. The following subsections describe the 
details for individual cases. 

3.1. Known binaries (Figs^ and^) 

Among seven known binaries with separations of 0.1-1.5 
arcsecs, displacement in Ha emission was detected in all 
of them but SR 20. Some of them also show displacement 
in other features including Fen 6456 A, Fen 6516 A, Hei 
6678 A, and Li I 6708 A. 

The position angles and the angular scales of the dis- 
placement are summarized in Table 2. All the observed 
displacements lie at the binary position angle, and this 
tendency can be attributed to the two binary components 
having different spectra (Bailey 1998b). The exception is 
Ha emission in S CrA and AS 205, in which the direction 
of the displacement differs from the binary position angle 
beyond the l-cr level. A likely explanation for S CrA is that 
the two binary components were not equally within the 
slit due to their wide binary separation. Indeed, infrared 
speckle observations give a binary separation of 1.4 arcsecs 
in this object, which is comparable to the seeing size for 
our observations (see Table 1). One might think that the 
observed difference of the position angle is due to orbital 
motion, although it is unlikely for the following reasons: 
(1) the binary position angle in S CrA increases with time 
(Bailey 1998b) while our results show a smaller position 
angle than previous measurements of the binary position 
angle by Ghez et al. (1997a), and (2) a central mass of 
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Table 1. Log of the observations 

Object HBC0 typ<|] V Distance| Date seeing Photons^] 1-ct Detection Limit 

(mag) (pc) (") (xlO 4 ) (mas) 



SZ Cha 


566 


K0 


12.0 


170 


1999. 


.7 


.2 


2.8 


2 


12.6 


Th 12 


605 


M0 


11.9 


190 


1999. 


.7 


.2 


2.1 


6 


2.9 


MWC 300 




B 


10.5 


1550 


1999. 


7. 


.2 


<1.9p| 
u 


4 


3.6 


T CrAfl 


290 


F0 


13.4 


130 


1999. 


.7 


.2 


<2.6 e 


4 


4.6 


S CrA 


286 


K6 


11.1 


130 


1999. 


.7 


.2 


<2.9 e 


2 


6.7 


LkHa 332-20 


244 


K2 


11.2 


170 


1999. 


.7 


.3 


1.2 


9 


1.3 


CS Cha 


569 


K5 


11.6 


170 


1999. 


7. 


.3 


1.2 


8 


1.6 


Sz 77 


603 


M0 


12.5 


190 


1999. 


7. 


.3 


0.9 


8 


1.5 


SR 20 


643 


<K7 


14.1 


140 


1999. 


7. 


.3 


1.4 


2 


3.9 


LkHa 118 


281 


B5 


11.1 




1999. 


.7 


.3 


1.8 


6 


3.3 


R CrA 


288 


A5 


10.7 


130 


1999. 


.7 


.3 


4.5 


2 


15.3 


LkHa 332-21 


247 


G8 


11.0 


170 


1999. 


7 


.4 


1.0 


14 


1.1 


Sz 41 


588 


K0 


11.6 


170 


1999. 


7 


.4 


1.1 


7 


1.9 


Sz 68 


248 


K2 


10.3 


190 


1999. 


7 


.4 


1.2 


23 


0.9 


GQ Lup 


250 


K7 


11.4 


190 


1999. 


7. 


.4 


1.3 


7 


2.1 


VV Ser 


282 


B-A 


11.9 


440 


1999. 


7 


.4 


1.2 


8 


2.6 


VV CrA 


291 


K7 


13.0 


130 


1999. 


7. 


.4 


1.2 


1 


5.0 


VW Cha 


575 


K2 


12.5 


170 


2000. 


.7. 


.21 


1.7 


1 


8.2 


AS 205 


254 


K5 


12.4 


150 


2000. 


7. 


.21 


2.0 


2 


5.6 


TW Hya 


568 


K7 


10.9 


56 


2000. 


7. 


.21 


1.3 


4 


2.2 


AS 353A 


292 




12.5 


300 


2000. 


.7 


.21 


2.3 


1 


7.7 


L CrA / 


290 


F0 


13.4 


130 


2000. 


7. 


.21 


<2.6 e 


2 


4.6 


V4046 Sgr 






10.5 




2000. 


7. 


.21 


2.8 


4 


4.5 


CU Cha 


246 


AO 


8.5 


170 


2000. 


7 


.21 


1.8 


1-2 


1.9 


AS 209 


270 


K5 


11.5 


150 


2000. 


7 


.21 


2.6 


3 


4.6 


HO Lup 


612 


Ml 


12.9 


190 


2000. 


7 


.23 


1.1 


0.5 


6.7 


CL Cha 


570 


K7 


12.4 


170 


2000. 


7. 


.23 


1.3 


2 


2.5 


BF Cha 








170 


2000. 


.7. 


.23 


1.3 


3 


2.8 


AS 206 


259 


K6-7 


12.9 


140 


2000. 


7 


.23 


1.4 


1 


5.8 



a the number in the Herbig& Bell (1988) catalogue 

b references — Herbig & Bell (1988), Leinert et al. (1997) 

c references — Bertout et al. (1999), de Geus & Burton (1991), de Zeeuw et al. (1999), Herbig 
& Jones (1983), Hillenbrand et al. (1992), Leinert et al. (1997), Marraco & Rydgren (1981), 
Whichmann et al. (1998) 

d for each three pixels (corresponding to 0.4 A) along the wavelength. 

e these objects show asymmetric profiles, suggesting that a stellar companion contributes to the 
measured FWHM 

* this object was observe twice to investigate the time variation of the spatial structure reported 
by Bailey (1998a) 



30-40 M© would be required to explain the difference be- 
tween the two position angles, not plausible as a mass for 
a PMS star (£8 M Q — Stahler & Walter 1993). On the 
other hand, the different position angles in AS 205 can be 
explained by the orbital motion of the system. Indeed, the 
difference between the two position angles suggest a mass 
of the primary star of 2.4 M© assuming a circular orbit, 
consistent with that derived from the stellar luminosity 
and spectral type (1.5-2.5 M© — Liu et al. 1996). 

The displacement of Ha emission/absorption in Sz 68, 
VW Cha, S CrA, and AS 205 show angular scales of 0.2- 
0.4 times as large as the binary separation, while it is 
much smaller in HO Lup and Sz 41 (0.01-0.02 as large 
as the binary separation). The latter results could be due 
to their large binary separations (^1.5 arcsec), in which 



the seeing of the brighter component does not include the 
flux from the fainter component. We thus conclude that 
angular scales for Ha emission in all the objects are consis- 
tent with those obtained in other pre-main sequence stars 
by Bailey (1998b): that is, 0.07-0.7. The displacement in 
most of the objects indicates the direction for one of the 
binary components, while that in the Ha profiles in S CrA 
and AS 205 lies in both directions: these profiles show cen- 
tral absorption and wing emission features, in which the 
position is displaced in the opposite directions. 

The angular scales of the displacement in other emis- 
sion/absorption lines are 0.02-0.05 per binary separation, 
much smaller than that in Ha emission. These small scales 
are attributed to smaller line-to-continuum ratios, which 
allow larger contribution of the continuum flux and reduce 



Takami et al.: Southern pre-main sequence stars 



5" 15 r 

J 5 h 

S jzf 

r WDO r 
Can 
I * n 

~100 r 
300 - 

E^lOO r 



VW Cha : Ho 8563A 




* t * 





0.3 p 

~ g L VW Cha : Hel B67BA 

e 1 f 
c 
" 50 F 




-50 

~ioo r 



i I i 



■J. o s-^yV^Af^fy^^^ 



-1000 



1000 



4V (km B _l ) 



1 '0 

I 6 

"308 
tna 

igsoo 

■"400 
200 

*? 



S CrA : Ha 6563A 




I I I 



i8 



*g B0 

wi o 

w 

-20 




i i i 



-1000 



1000 



aV (km B _I ) 



Fig. 1. Examples of intensity and position spectra for known binaries. Intensity spectra are normalized to the contin- 
uum. 



the positional displacement. In S CrA, the displacement 
indicates different directions between the lines: two Fen 
lines lie in the direction of the primary, while He I and 
[S il] lines lie in the direction of the secondary. 

Positional displacement of forbidden lines in PMS stars 
is often attributed to a jet/wind (e.g., Hirth, Mundt, & 
Solf 1997; Takami et al. 2001,2002), however, this expla- 
nation is not likely for the [Sll] lines in S CrA. Indeed, 
optical polarization of T Tauri stars is very often perpen- 
dicular to jets within a few degrees (Menard & Bastien 
1992), although that in S CrA was measured to be 151— 
156 degree (Bastien 1985), nearly parallel to that of the 
positional displacement in the [S n] emission. We thus con- 
clude that the displacement of the [S n] lines of S CrA is 
due to a binary companion. 

3.2. T CrA (Fig. |j 

Infrared speckle observations by Ghez et al. (1997a) did 
not show the presence of a stellar companion, whereas 
the position spectra show displacement similar to that in 
known binaries, as first reported by Bailey (1998b). His 
and our results show angular scales for the displacement of 
76±5, 140±9, 122±5, and 113±3 milliarcsec in 1996, 1997, 
1999, and 2000, respectively, indicating a lower limit for 
the binary separation of 140 milliarcsec. Such a separation 



is larger than the detection limit of infrared speckle, thus 
the non-detection of the companion by Ghez et al. (1997a) 
implies that the companion has a if -magnitude fainter 
than 10.5. 

Fig. ^ shows the position angles of the displacement 
measured in our four observing runs. The angle increases 
with time, corresponding to ^2 degree/yr, consistent with 
the orbital motion of the stellar companion. From this 
result, we estimate a binary separation of ~260 milliarcsec 
adopting a primary stellar mass of 0.8 M (Hillenbrand et 
al. 1992) and assuming a circular orbit. In this case, the 
angular scales observed in the four runs would correspond 
to 0.29-0.54 per binary separation, consistent with those 
observed in known binaries (§3.1). On the other hand, time 
variation of the observed angular scale is not systematic 
as expected for binary motion, and such results could be 
attributed to variability of the emission feature between 
binary components (Bailey 1998b). Indeed, Ha emission 
observed by Bailey (1998b) and this work exhibit different 
shapes of profiles, supporting this explanation. 

3.3. MWC 300 (Fig. @J 

Positional displacement was detected in the central ab- 
sorption feature of the Ha emission, with an angular scale 
and a position angle of 17±4 milliarcsec and ^165°, re- 
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Fig. 2. Intensity profiles and positional displacement in emission/absorption features of known binaries VW Cha, S 
CrA, Sz 68, and AS 205. Small crosses in each map correspond to the velocity components shown in the intensity 
profile. A large cross in each figure indicates the centroidal position of the continuum source. A dot-dashed line indicates 
the direction of the fainter companion in the it"- band continuum (Ghez et al. 1993,1997a). The intensity profiles are 
normalized to the continuum. 
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in jets/winds (e.g., Takami et al. 2001,2002), where the 
displacement is observed at negative or both velocities. 

Previous infrared speckle observations failed to detect 
this companion (Leinert, Richichi, & Haas 1997; Pirzkal, 
Spillar, & Dyck 1997). Corporon & Lagrange (1999) ar- 
gue that MWC 300 is a spectroscopic binary, although 
our results cannot be explained as such: a combination of 
a lower limit for the binary separation provided by our 
results, the distance to the star (15.5 kpc — Leinert et al. 
1997), and a typical mass of Herbig AeBe stars (5 M Q ) 
suggests an orbital period of more than 2 x 10 3 yrs, much 
larger than the duration of the observations by Corporon 
& Lagrange (1999). 



Fig. 4. Position angles of the displacement measured in T 
CrA. A dot-dashed line shows the regression line with an 
inclination of 1.65° yr _ . 

spectively. Similar positional displacement was also ob- 
served in the known binary S CrA; thus, it is likely that 
the observed displacement in MWC 300 is due to a stellar 
companion. One might think that the observed displace- 
ment could be explained by a jet/wind, however, it is not 
likely since the displacement in this object is observed 
only at positive velocities, different from those expected 



3.4. CS Cha (Figs @ and @J 

The displacement in Ha emission shows a bipolar struc- 
ture: the blueshifted component extends towards the 
north-east while the redshifted component extends to- 
wards the south-west. The direction of the displacement 
is nearly perpendicular to that of the optical polarization 
(139° — Drissen, Bastien, & St. -Louis 1989), which is very 
often perpendicular to that of jets (Menard & Bastien 
1992, and references therein). Thus, it is likely that the 
bipolar structure in CS Cha is due to outflowing gas, as 
Takami ct al. (2001) argue for RU Lup. 
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Table 2. Positional displacement in known binaries 



object 


separation^ 


P.A." 






line 




(mas) 


(deg.) 








Sz 68 


107 ± 7 


245 ± 


7 


Hq 


6563 A 


Sz 41 


1500 ± 800 


150 ± 20 


Hq 


6563 A 


HO Lup 


1490 ± 70 


DO ± 


n 

z 


Hq 


6563 A 


VW Cha 


660 ± 30 


184 ± 


2 


Hq 


6563 A 










He I 


6678 A 


S CrA 


1410 ± 60 


157 ± 


2 


Hq 


6563 A 










Fen 


6456 A 










Fen 


6516 A 










He i 


6678 A 










[Sn] 


6716 A 










[Sill 


6731 A 


AS 205 


1317 ± 2 


211.87 ±0.( 


)8 


Hq 


6563 A 










Fen 


6456 A 










Fen 


6516 A 










He i 


6678 A 










Li i 


6708 A 



a observed by Ghez et al. (1995,1997a) 



observed angular scale observed P.A. 

(milliarcsec) (per binary separation) (deg.) 



21 ± 1 


0.20 


248.2 ±0.9 


22 ± 3 


0.01 


186 ±3 


24 ± 9 


0.02 


213 ± 1 


250 ± 12 


0.38 


2.5 ±0.5 


93 ± 12 


0.14 


4±3 


513 ± 9 


0.36 


153.1 ±0.3 


111 ± 9 


0.08 


331.1 ±0.5 


47 ± 9 


0.03 


323 ±3 


49 ± 9 


0.03 


338 ±4 


57 ± 9 


0.04 


150 ±5 


31 ± 9 


0.02 


146 ± 2 


57 ± 9 


0.04 


148 ±3 


239 ± 8 


0.18 


33.5 ±0.4 


67 ± 8 


0.05 


202 ±4 


61 ± 8 


0.05 


36 ±3 


60 ± 8 


0.05 


37 ±2 


34 ± 8 


0.03 


34 ±4 


81 ± 8 


0.06 


28 ±3 



To determine the true spatial scale of the bipolar flow, 
we removed contamination by the continuum using the 
line-to-continuum ratio at each wavelength (see Takami 
et al. 2001 for detail), and plot the results in Fig. ||. 
The figure shows an angular scale of 8-10 milliarcsec, and 
asymmetric velocity field between the blueshifted and red- 
shifted component: the displacement is observed at veloc- 
ities of —100 to —20 km s _1 in the blueshifted component 
and ±100 to ±250 km s _1 in the redshifted component. 
The displacement along the outflow axis does not show 
clear positive correlation with the velocity as observed 
in RU Lup by Takami et al. (2001). Thus, the velocity 
field observed in this object cannot be simply attributed 
to the acceleration of the outflow, in which the displace- 
ment should increase with the velocity. Our result could 
be due to a contribution of magnetospheric mass accre- 
tion columns (see, e.g., Muzerolle, Calvet, & Hartmann 
1998,2001) or an inner hot wind (Beristain, Edwards, & 
Kwan 2001; Alencar & Batalha 2002) to the observed Ha 
flux. 

Forbidden lines in PMS stars often exhibit only 
blueshifted components due to obscuration of the counter 
flow by a circumstellar disc. (Eisloffel et al. 2000, and refer- 
ences therein). On the other hand, Ha outflows in CS Cha 
have both blueshifted and redshifted components. This re- 
sult can be explained if a gap or hole in the circumstellar 
disc allows the redshifted components to be seen (Takami 
et al. 2001). Indeed, the infrared spectral energy distribu- 
tion (SED) in CS Cha shows no substantial excess emis- 
sion less than 20 /im (Gauvin & Strom 1992), indicating 
a lack of dust with temperatures of more than 200 K. If 
we adopt a stellar luminosity of 2.7 L Q (Gauvin & Strom 
1992), this temperature corresponds to a disc radius of 



6 AU, consistent with the spatial scale of the observed 
redshifted Ha outflow, i.e., ~1.4 AU. 

3.5. AS 353 A (Fig. gj 

A positional displacement was observed within the 
blueshifted absorption of the Ha profile, at a velocity 
range between —200 and —300 km s _1 . This velocity co- 
incides with that of the high-velocity component in the 
forbidden line profiles, peaking at —300 km s _1 (Hamann 
1994) and considered to arise from a jet close to the star 
(Eisloffel et al. 2000 for review). In addition, the position 
angle of the displacement of 111±18° is coincident with 
that of a HH object associated with this star, HH 32C 
(103° — Solf, Bohm & Raga 1986). We thus conclude 
that a jet close to the star is responsible for the observed 
displacement. One might think that a binary companion 
could also explain the observed displacement, however, it 
is not likely since the emission component in the Ha pro- 
file does not show any displacement as observed in known 
binaries (see §3.1). 

Fig. H shows that the displacement has an angular scale 
of 70 milliarcsec, corresponding to ^20 AU at the star. 
The inclination of the extended jet is measured to be 70° 
from the plane of the sky (Curiel et al. 1997); thus, the cen- 
troidal position of the emission line region is presumably 
located more than 60 AU apart from the star. Position 
spectra in this object do not show a displacement corre- 
sponding to the counterflow, as observed in RU Lup and 
CS Cha. Such results may be due to obscuration of the 
flow by a circumstellar disc, as observed in forbidden lines 
in many T Tauri stars. 
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Fig. 3. Intensity profile and position displacement of the Ha emission in T CrA, MWC 300, CS Cha, AS 353A, and 
R CrA. The intensity profiles are normalized to the continuum. The results of R CrA in 1996 was obtained by Bailey 
(1998b). 



3.6. R CrA (Fig. gj 

The positional displacement observed in Ha emission dif- 
fers from any other objects described above. In 1996, the 
bluer component at —210 to 40 km s _1 is displaced with 
an angular scale of 120 milliarcsec and at a position angle 
of ~ 220°, while the redder component at 110 to 600 km 
s _1 docs not show clear displacement. In 1999, the posi- 
tion of these components are displaced at position angles 
of -200° and -130° with angular scales of 100 and 70 
milliarcsecs, respectively. Neither a stellar companion or 
a monopolar/bipolar outflow can explain these results. 

Such results could be explained by a combination of 
two companions, or that of a stellar companion and an 
outflow. However, the former explanation is not likely for 
the following reason. Recent surveys show that young mul- 
tiple systems are invariably hierarchical: i.e., with differ- 



ent orbital scales by a factor of at least 3 and typically 
more than 10 (e.g., Simon et al. 1995; Ghez et al. 1997a; 
Kohler et al 2000). This implies that two stellar com- 
panions cannot simply explain our results in R CrA, in 
which displacement in two directions show similar angu- 
lar scales. A combination of a stellar companion and an 
outflow is thus more likely to explain the observed dis- 
placement. Emission lines from outflows often exhibit only 
blueshifted components (e.g., Hirth, Mundt, & Solf 1997): 
thus, a likely explanation is that the displacement in the 
bluer emission is due to an outflow, while that in the red- 
der emission is due to a stellar companion. 

Hillenbrand et al. (1992) estimate a primary stellar 
mass of 0.8 M Q in this object, and our results provide 
a lower limit for the binary separation of 8 AU. These 
parameters provide the timescale of the orbital motion of 
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Fig. 5. (left) Spatial distribution of positional displacement in Ha emission of CS Cha. Contamination by the contin- 
uum was removed using the line-to-continuum ratio at each wavelength: thus, the plot shows the true angular scale 
of the Ha emission. Open and filled circles show blue- and rcdshiftcd components, respectively, and each circle in the 
map corresponds to that in the intensity profile shown below. A large cross and a dashed bar at the centre show the 
centroidal position of the continuum and the direction of optical polarization, respectively. The scale of the profile 
is normalized to the continuum flux, (right) Correlation between velocity and radial distance of the displacement. A 
large cross shows the centroidal position of the continuum. 



more than 24 yr, much larger than the duration of our 
observations. This fact suggests that the observed time 
variation in the redshifted emission is due to the variation 
of the spectra in two binary components, as observed in 
T CrA (§3.2), rather than binary motion. 

3.7. VV CrA (Fig.® 

Positional displacement was observed in Ha and several 
forbidden lines, and these are attributed to the following 
three flows: (1) flow A, which appears in the [S n] and [O i] 
emission at velocities of —100 to 20 km s _1 ; and (2) flow B, 
which appears in the [N n] and [S n] emission at velocities 
of —270 to —100 km s _1 ; (3) flow C, which appears in the 
[Ol] and Ha emission at velocities of —350 to —200 km 
s . The position in the flow monotonically changes with 
velocity in the [Sn] emission of flow A and B, and that in 
the [Ol] line in flow C. 

The displacement in the lowest velocities of flow A is 
associated with the primary star shown in Fig. ^, indicat- 
ing that the flow is driven by the star at the centroidal con- 
tinuum position. The [S n] profiles exhibit distinct compo- 
nents at high and low velocities, as observed in many other 
stars, and the latter component is located within 0.2 arc- 
sec of the star, corresponding to spatial scales of < 25 AU. 
This spatial scale is similar to those observed in other T 
Tauri stars (e.g., Hirth, Mundt, & Solf 1997; Takami et 
al. 2001,2002), supporting the idea that the flow is driven 
by this PMS star. On the other hand, in flow B, the dis- 
placement at the lowest velocities lies in the north-east, in 
the direction of a companion with a separation of 2 arc- 
sec (Graham 1992; Reipurth & Zinnecker 1993; Ghez et 



al. 1997a). Thus, it is likely that flow B is driven by this 
companion, known as an "infrared companion" (Koresko, 
Herbst & Leinert 1997). The [Nil] emission in this flow 
is located further from this star than the [SII] lines, as 
observed in small-scale jets in many other T Tauri stars 
(Hirth, Mundt, & Solf 1997), and such results support this 
conclusion. 

Previous infrared observations by Graham (1992) and 
Ghez et al. (1997a) do not show the presence of a young 
star associated with flow C. In this flow, the [O i] line has a 
width of a few hundred km s _1 , suggesting the coexistence 
of shocks with different velocities. In addition, the line 
with higher excitation energy, Ha, arises from the highest 
velocities. These results can be attributed to a bow shock 
in which the Ha and [O i] emission arise from its apex and 
wakes, respectively. The extension of the [O i] emission 
suggest that the driving source lies in the west, in which a 
few IRAS sources were identified by Wilking et al. (1992), 
although it is located far (>2 arcmin) from VV CrA. Thus, 
the driving source of flow C is not clear. 

4. Ha profiles in binary components 

We can derive the intensity spectrum for each component 
of known binaries from the intensity and position spectra, 
the binary separation measured by infrared speckle tech- 
niques, and an assumed continuum flux ratio (see Bailey 
1998a; Garcia, Thiebaut, & Bacon 1999). The centroidal 
position of the seeing in a PMS binary is described as 
follows: 

_ ^A(prim) ^ A (sec) ^ 

X A — j x (prim) i" j ■ x -(sec)i 
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Fig. 6. Positional displacement in VV CrA. Contamination by the continuum is removed as carried out for Fig. p[ 
Small open squares, filled and open circles, and stars correspond to velocity components in the profiles of [Nil] 6583 A, 
[Sn] 6731 and 6736 A, and Ha, respectively, at the bottom-left of the figure. Small open triangles and crosses are those 

for [Oi] 6300 A but different velocities 330 to —200 km s _1 and —180 to 10 km s _1 , respectively. Different grey 

tones indicate the displacement in flow A-C, respectively (see the text for detail). A large cross shows the centroidal 
position of the continuum. The dot-dashed arrows show the variation of the velocities in the flows from low to high 
values. The displacement in [O i] 6363 A is not plotted, since it is consistent with that in [O i] 6300 A, which has the 
same upper transition level. The intensity profiles in the bottom-left of the figure are normalized to the continuum. 



or 

Ad\ = (2) 



where is the centroidal position of the seeing, and 
x (prim) an d X( sec ) are positions of the primary and the 
secondary stars, respectively; I\, ^(prim); ^a(scc) are inten- 
sities of the whole system, the primary star, and the sec- 
ondary star, respectively; Ad\ is the separation between 
the primary star and the centroidal position of the seeing; 
d is the binary separation projected on the sky. Positional 
displacement in an emission/ absorption feature is mea- 
sured relative to the centroidal position of the continuum 



emission, thus its spatial scale is described as follows based 
on equation (2): 

Ad' x = Ad\ - A4ont = ( ^1 _ htp&L) d) (3) 

\ Jcont / 

where Ad' x is the observed displacement in a position 
spectrum; Ad cont , I C ont and / con t( se c) are the distance to 
the centroidal continuum position from the primary star, 
the total luminosity of the continuum, and that of the 
secondary star, respectively. From this equation, the flux 
from each star is obtained as follows: 

/a(scc) - vr^r + —) h - [ttj + —) /a ' 

^A(prim) = I\-I\(sec), (4) 
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where / is the continuum luminosity of the secondary star 
relative to the primary. 

Fig. ^ shows the profiles of the Ha emission obtained 
from equation (4) for six of our objects with known sep- 
aration. Despite the uncertainty of the continuum flux 
ratio, the figure shows that the two binary components 
have strikingly similar profiles to each other — similar 
lincwidths and the absorption feature at similar veloci- 
ties. The figure also shows that the line-to-continuum ra- 
tio is higher at the secondary than the primary in all the 
objects. Although the profiles in HO Lup do not clearly 
show this tendency, the observed positional displacement 
lies in the direction of the fainter component (Reipurth 
& Zinnecker 1993), agreeing with this conclusion. These 
results also imply that the position of the Ha feature in 
all the objects is primarily displaced in the direction of 
the component fainter at optical wavelengths. 

5. Discussion 

5.1. Detectability of close (<0.1 arcsec) binaries using 
spectro-astrometry 

Recent high-resolution techniques, such as infrared speckle 
observations and the Hubble Space Telescope, have re- 
vealed a difference in binary populations between star- 
forming regions. These observations show that the binary 
frequencies in low-mass star-forming regions are higher 
than that of solar-type field stars with the same range of 
binary separation (e.g., Ghez et al. 1993; Leinert et al. 
1993; Ghez et al. 1997a). Alternatively, the frequency ap- 
pears to be similar to that of the field stars in a massive 
star-forming region, the Orion Trapezium cluster (Prosser 
et al. 1994; Padgett et al. 1997; Petr et al. 1998). Such 
a different binary population between regions could be 
caused by initial conditions, such as stellar density or 
cloud temperature (see Mathieu et al. 2000 for review). 

However, the traditional observing techniques have not 
been able to cover the complete range of binary separa- 
tions in various star-forming regions. Direct imaging and 
infrared speckle have allowed detection of PMS binaries 
down to separations of ^100 milliarcsecond, correspond- 
ing to ~15 AU with periods of 50 years or more in the 
nearest star-forming regions. On the other hand, spectro- 
scopic measurement has been capable of detecting binaries 
with periods of a few days and separations of well under an 
AU. These techniques cannot cover the intermediate sepa- 
rations (0.1-10 AU), leading to significant incompleteness. 
Lunar occultations have explored this gap in some cases, 
however, this technique can measure the separation only in 
one-dimension in the sky, and applies for a limited num- 
ber of star-forming regions as well as relying on chance 
occultations. 

Spectro-astrometric observations by Bailey (1998b) 
and this work suggest that this technique is capable of 
discovering PMS binaries with these intermediate separa- 
tions. These observations detect 10 out of 12 known bina- 
ries with separations of 0.07-1.5 arcsec, and a few new bi- 



nary candidates which have not been detected using other 
techniques. The observed angular scales in known binaries 
(0.07-0.7 times as large as the binary separation — §4) 
and positional accuracy achieved (down to 1 milliarcsec) 
suggest that the technique should be able to detect bina- 
ries with separations down to ~10 milliarcsec at the 3-cr 
level. This detection limit is comparable to that of lunar 
occultations (see Richichi et al. 1994; Simon et al. 1995), 
although spectro-astrometry can easily achieve better lim- 
its by detecting a larger number of photons (see equation 
1), and applies for objects with any coordinate and binary 
position angles. 

Equation (3) shows that the detectability of binaries 
using this technique depends on the line and continuum 
flux ratios between binary components. In addition, the 
position spectra in Ha emission show different shapes from 
the intensity profiles (see Fig. 1), indicating that the de- 
tectability also depends on a slight difference in the line 
profiles for the primary and secondary stars. These facts 
may be susceptible to more complicated selection effects 
than those of direct imaging, infrared speckle, and lunar 
occultations. To investigate this effect, we plot the rela- 
tion between the angular scale of the displacement, Ha 
equivalent width of the whole system, and relative K mag- 
nitude between binary components in Fig. ||. The figure 
shows that the displacement has no correlation with the 
Ha equivalent width, while there is a clear correlation with 
the relative K magnitude. These results suggest that our 
technique has a selection effect similar to those of infrared 
speckle (Ghez et al. 1993) and lunar occultations (Simon 
et al. 1995); i.e., dependent on the relative K magnitude 
and independent of the line flux. The latter correlation 
between the displacement and relative K magnitude can 
be explained as follows: if the secondary star has a bright 
stellar / disc luminosity, it is likely to have bright Ha emis- 
sion (see Cabrit et al 1990) and bright K band emission, 
thereby providing a small difference in the K magnitude 
from that of the primary and a large displacement in Ha 
emission towards the secondary. 

From the regression line in Fig. ||, we estimate the 
detection limit of binary separation using our technique as 
a function of the relative K magnitude (Fig. ^). The same 
figure also shows binary parameters discovered by lunar 
occultations with separations of less than 0" .1 (Richichi et 
al. 1994,1996,1997; Simon et al. 1995). The figure indicates 
that we could detect 90 % of binaries discovered by lunar 
occultations, and our detection limit appears to be slightly 
better than that of the lunar occultation method for both 
binary separation and relative K magnitude. Such a high 
sensitivity for relative K magnitude is also followed by 
a discovery of a companion candidate in T CrA, which 
Ghez et al. (1997a) failed to detect despite a sufficient 
separation (>0.14 arcsec — see §3.2). 

In two separated studies, we failed to detect the com- 
panions in T Tau (Bailey 1998b) and SR 20 (this work), 
which were detected by Ghez et al. (1993) using infrared 
speckle. The non-detection for the former object could be 
due to the fact that the secondary star is an "infrared com- 
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Fig. 7. Profiles of the Ha emission for each binary component of known binaries Sz 68, Sz 41, HO Lup, AS 205, VW 
Cha, and S CrA. Solid, dot, and dot-dashed profiles were obtained by assuming the flux ratios of /=0.2, 0.4, and 
0.8, respectively, at the continuum level. In this figure, the "primary" and "secondary" components arc defined based 
on the continuum flux at 6500-6800 A, different from Fig. || and Table 2. The intensity scale is normalized to the 
continuum for each binary component. 



panion" and faint at optical wavelengths (Koresko et al. 
1997). Such companions form only ~10% (White & Ghez 
2001) or less (Ghez et al. 1997a; Kohler et al. 2000) of 
the entire PMS population. Another possible explanation 
for non-detection is that Ha emission is faint or absent at 
the companion of T Tau and SR 20, thereby not allowing 
the positional displacement towards the secondary. Again, 
such binary systems may not significantly affect the en- 
tire population for PMS binaries: recent observations have 
shown that companions of classical T Tauri stars (CTTSs) 
are usually CTTSs, and a combination of a CTTS and a 
weak-line T Tauri star (WTTS) is rare (^20% — Prato 
& Simon 1997; Duchene et al. 1999). 

All the results described above suggest that we can 
improve the previous statistics of PMS binary popula- 
tions using spectro-astrometry. Note that most of the 
known binaries in our sample have bright Ha emission 
(£ , M / H Q>10 A), implying that the detection of PMS bi- 
naries using our technique could be efficient only for such 
objects as CTTSs. Indeed, equation (3) suggests that the 
displacement could be shorter for binaries with small line- 
to-continuum ratios. Although Fig. || does not show such 
a tendency, it could appear if we included the objects with 
fainter Ha emission. This fact implies the detectability of 



binaries using our technique could be worse for objects 
such as WTTSs. 

5.2. Binary frequency in our sample of the nearest 
southern star-forming regions 

The observed targets shown in Table 1 include 22 PMS 
stars in nearby (130-200 pc) star-forming regions: includ- 
ing 7 known binaries with separations of 0.07-1.5 arcsec 
(§3.1), and two "spectro-astrometric" binaries in T CrA 
(§3.2) and R CrA (§3.6). Among these binaries, R CrA 
may have a separation of less than 0.1 arcsec, correspond- 
ing to ^15 AU. In contrast, lunar occultations surveys by 
Simon et al. (1995) detected 8 binaries out of 47 PMS 
stars in Taurus and 5 out of 35 PMS stars in Ophiuchus 
at these short (5^15 AU) separations. Is the binary fre- 
quency at these separations smaller in our sample than 
the Taurus and Ophiuchus star-forming regions? 

To investigate this, we determine the binary frequency 
of our sample. We exclude Sz Cha from this statistical 
study due to a remarkably low detection limit (13 mas). 
The objects in which a companion has not been detected 
have a typical uncertainty in their measured displacement 
of 3 milliarcsec. Given that the typical observed angular 
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Fig. 8. Relation between the angular scale of the displace- 
ment and Ha equivalent width of the whole system (up- 
per), and relative K magnitude between binary compo- 
nents (lower). The open and filled circles show the results 
for Herbig Ae/Be stars and classical T Tauri stars, re- 
spectively. The regression line (dashed) is shown in the 
lower figure. Data for the displacement and Ha equivalent 
widths are obtained from Bailey (1998b) and this work, 
while the relative K magnitudes are obtained from pre- 
vious infrared speckle observations by Ghez et al. (1993), 
Ghez et al. (1997a), and Leinert et al. (1997). 
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Fig. 9. The detection limit of PMS binaries using spectro- 
astrometry. The hatched area shows the parameter space 
in which we cannot detect binaries exceeding the 3-<r level. 
This limit is estimated assuming a positional uncertainty 
of 1 mas in the position spectra. The filled and open cir- 
cles, and crosses show the binaries discovered by lunar 
occultations by Simon et al. (1995), Richichi et al. (1994), 
and Richichi et al. (1996, 1997), respectively. 



Table 3. Binary population in star-forming regions 
- separation = 15-300 AUf] - 



region 


binary fraction 


reference 




la/urus 


n no 
u.o t ztu.uy 


Grhez et al. 


l^ryyo ) 




0.32±0.08 


Simon et al. 


(1995) 


Ophiuchus 


0.14±0.06 


Simon et al. 


(1995) 


Oph-Sco 


0.29±0.12 


Ghez et al. 


(1993) 


Sco-Cen 


0.26±0.05 


Kohler et al. 


(2000) 


Cha/Lup/CrA 


0.31±0.08 


Ghez et al. 


(1997a) 


field stars 


0.16±0.04f] 


DM9l|] 




mir sflmnlp in 

U HI odlilLJirj III 








southern SFRs 


0.41±0.16 


this work 







separation = 5-15 


Air - 




region 


binary fraction 


reference 




Taurus 


0.06±0.04 


Simon et al. 


(1995) 


Ophiuchus 


0.06±0.04 


Simon et al. 


(1995) 


field stars 


0.06±0.02 b 


DM91 C 




our sample in 








southern SFRs 


0.05±0.05 


this work 





a assuming the distance of all the star-forming regions 
of 150 pc for simplicity 

b adopting the factor by Leinert et al. (1993) to convert 
the orbital period to the separation 

c Duquennoy & Mayor (1991) 



scale per binary separation is 0.3 for Ha emission (see §3.1 
and Table 1), this implies that a companion has not been 
detected down to 30 milliarcsec (3er), equivalent with ^5 
AU in the observed star-forming regions. Table 3 shows 
binary populations in our sample with two different sep- 
aration ranges (15-300 AU and 5-15 AU), together with 
those in various star-forming regions and field stars. An 
accurate binary separation is not known for R CrA, and 
this fact is included in the uncertainty of the tabulated 
values. The table shows that our sample has a larger bi- 
nary frequency for separations of 15-300 AU than that in 
field stars. Such a result is consistent with those derived in 
other low-mass star-forming regions. On the other hand, 
all the frequencies tabulated for the shorter separations 
have large statistical uncertainties, implying that the pro- 
vided separation range is insufficient to investigate the 
difference of binary frequencies between star-forming re- 
gions. Observations of 60 stars with an uncertainty of 1 
mas in their position spectra could allow us to investigate 
it in detail: we could detect 8±3 companions with these 
separation ranges if the binaries in the observed regions 
had the same number distribution of separations as that 
in solar-type main-sequence stars observed by Duquennoy 
& Mayor (1991). 
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5.3. Implications for the binary formation mechanism 



Table 4. Summary of Ha outflows 



Approximately two-thirds of all solar-type main-sequence 
stars are members of multiple systems (Duquennoy & 
Mayor 1991), and it may be even higher in low-mass star- 
forming regions (see Mathieu et al. 2000). Despite the ap- 
parent ubiquity of binaries, astronomers have not yet iden- 
tified the primary mechanism by which they form. Possible 
formation mechanisms include (see Hartigan, Strom, & 
Strom 1994a; Ghez, White, & Simon 1997b): (1) frag- 
mentation of a collapsing molecular core, (2) gravita- 
tional instability in an equilibrium disc, and (3) dynam- 
ical capture. The secondary mass distribution in field 
stars is remarkably similar to the mass function of pri- 
mary stars with low-to-high masses (Fischer & Marcy 
1992; Duquennoy & Mayor 1991; Abt, Gomez, & Levy 
1990), supporting protostellar "capture" as the mecha- 
nism by which binaries form. On the other hand, stellar 
ages, infrared excess, and emission line features in PMS 
binaries indicate that these are at the same evolutionary 
stage, supporting the scenario of core fragmentation (e.g., 
Hartigan et al. 1994; Ghez et al. 1997b; Prato & Simon 
1997; Duchene et al. 1999; White & Ghez 2001). 

Our results in §4 show that the two binary compo- 
nents have strikingly similar Ha profiles to each other. 
Such results contrast to a rich variety of profiles among 
PMS stars (see Reipurth, Pedrosa, & Lago 1996). The va- 
riety of these profiles is due to interplay between the accre- 
tion flow and a wind closest to the source: thus, our results 
suggest that the two binary components tend to have simi- 
lar flow geometry and velocity for mass accretion/ejection. 
Such results could exclude the scenario of dynamical cap- 
ture and disc instability for binary formation, since (1) 
capture would provide a rich variety of profiles even be- 
tween binary components, and (2) disc instability could 
provide different evolutionary stages between binary com- 
ponents, thereby providing different inflowing/outflowing 
activities. Our results thus support the scenario of core 
fragmentation for binary formation, as well as other stud- 
ies for PMS binaries 

Accretion flows and/or winds in YSOs are likely to 
have non-spherical geometry (Shu et al. 2000; Konigl et 
al. 2000; Najita et al. 2000): thus, the similarity of the 
Ha profiles observed in binary components would also re- 
quire the condition that the axes of the inflows/outflows 
are aligned with each other. Such a condition is consistent 
with a previous polarimetric study by Monin, Menard, & 
Duchene (1998), who show that the position angles of po- 
larization are preferentially parallel in PMS binaries. On 
the other hand, Hale (1994) shows that the orbital and 
equatorial planes appear to be randomly distributed in 
main-sequence binaries with separations more than 30-40 
AU. Such a difference between PMS and main-sequence bi- 
naries could result from different formation mechanisms, 
or different initial conditions of the protostellar core. In 
the latter case, the results in PMS binaries could come 
from the pre-collapse material all rotating around a single 
axis, while those in main-sequence binaries could origi- 



object 

RU Lup 
CS Cha 



projected scale 
(AU) 

~1.5 



Z CMa > 7C(j 

AS 353A > 20 



velocity" 
(km s" 1 ) 

-200<u< -50 
200<u< 300 

-100<«< -20 
100<«< 250 

-700<v< -50 

-300<u<-200 



EW H! E 

(A) 

29.7 
28.3 

-B 



a sum for the blueshifted and redshifted emission in 
which positional displacement is observed 

b larger than that in Takami et al. (2001) due to the 
updated distance (Wichmann et al. 1998) 

c adopting the distance of 1150 pc provided by Claria 
(1974) 

d equivalent widths cannot be derived since the dis- 
placement is overlapped with the absorption feature 



nate from an the initial cloud core in which the angular 
momentum distribution has strong spatial variations (see 
the discussion of Bate et al. 2000). 

5.4. Physical Conditions of Jets down to AU Scales 

The results for the Ha outflows obtained by Bailey 
(1998b), Takami et al. (2001) and this work are summa- 
rized in Table 4. Among these objects, three of them are 
known as active PMS stars: Z CMa and AS 353 A are 
associated with extended jets (e.g., Poetzel, Mundt, & 
Ray 1989; Solf, Bohm, & Raga 1986) and RU Lup shows 
extremely high accretion luminosity (3-6 Lq — Lamzin 
et al. 1996). Such high activities have not been explic- 
itly reported for CS Cha, although the object shows high 
variability in Ha emission:- measured equivalent widths 
range from 13 A (Herbig & Bell 1988) to 54 A (Reipurth, 
Pedrosa, & Lago 1996) and 60 A (this work), and the ob- 
served line profile in Fig. |B| show a different shape from 
that observed by Reipurth et al. (1996). 

The velocities of blueshifted Ha outflows in RU Lup, 
Z CMa, AS 353A coincide with those of the high-velocity 
components in forbidden lines: —100 to —250 km s _1 in 
RU Lup (Takami et al. 2001); -600 to -800 km s" 1 in 
Z CMa (Corcoran & Ray 1997); -200 to -350 km s" 1 in 
AS 353A (Hamann 1994). This forbidden line emission is 
considered to arise from jets close to the star (see Eisoffel 
et al. 2000, and references therein): thus, this coincidence 
of velocities suggests that the Ha emission in these objects 
also arises from jets close to the source. Although the Ha 
emission in PMS stars are often considered to arise from 
the magnetospheric accretion columns, these results are 
not surprising for the following reasons. First, the ideal- 
ized magnetospheric accretion models by Muzerolle et al. 
(1998,2001) failed to reproduce permitted line profiles in 
active PMS stars, suggesting that the presence of another 
contributor to the emission. In addition, dozens of YSOs 
exhibit thermal radio emission, which is likely due to the 
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ionized gas in jets close to the source (see Anglada 1996). 
From the Ha luminosities in RU Lup and CS Cha, we can 
derive upper limits of the thermal radio luminosities of 
0.1-0.4 mJy kpc 2 at 3.6 cm (see Appendix A), consistent 
with those observed in other Class II objects (sec Bcltran 
et al. 2001)[]. 

Ha emission-line regions in Z CMa and AS 353A show 
spatial scales more than 50 AU (see Table 4 and §3.5), and 
there is growing evidence that the jets at these scales are 
heated by shocks in the internal working surfaces. Recent 
high-resolution imaging of PMS stars exhibit knotty, bub- 
bly, and bow-shape structures in jets at these scales (e.g., 
Kepner et al. 1993; Bacciotti et al. 2000; Dougados et al. 
2000). Lavalley-Forquet et al. (2000) carried out a detailed 
analysis of line ratios in DG Tau, and reach the conclu- 
sion that the jet is heated by shock, and not ambipolar 
diffusion nor turbulent viscosity in the mixing layer. Their 
conclusion is supported by Takami et al. (2002), who dis- 
cover He I 1.083 /im emission in the same jet: only shocks 
could explain the presence of such a high-excitation (~20 
eV) line. 

On the other hand, the same heating mechanism can- 
not simply explain the following results in RU Lup and 
CS Cha. Ha emission in these objects show spatial scales 
of 1-5 AU, much smaller than the distance between in- 
ternal working surfaces observed in PMS jets (^SlOO AU 
- see, e.g., Ray et al. 1996; Dougados et al. 2000). In 
addition, their line profiles do not exhibit a peak at the 
velocity of the jets, and indicate that Ha emission in these 
jets is continuously distributed from the accelerating re- 
gion. These results suggest the presence of another heat- 
ing mechanism in the inner (<10 AU) region, which could 
be (1) X-ray heating (Shang et al. 2002); (2) MHD-wave 
(e.g., Hartmann, Edwards, & Avrett 1982); (3) ambipolar 
diffusion (Safier 1993; Garcia et al. 2001); (4) turbulent 
dissipation in a viscous mixing layer (Binette et al. 1999); 
or (5) UV photoionization. Among these mechanisms, am- 
bipolar diffusion and viscous mixing layer models provide 
a lower ionization fraction than shocks, as indicated by the 
predicted [Nil]/ [OI] line ratios (see Lavalley-Fouquet et al. 
2000). Such a prediction is inconsistent with the spatial 
distributions of Ha and [OI] emission in RU Lup, which 
indicates a higher ionization fraction in the inner region 
than for the outer region heated by shocks. PMS stars 
often exhibit UV continuum excess as a result of mass 
accretion (see e.g., Calvet, Hartmann, & Strom 2000), al- 
though this could also be rejected for the following reason. 
The hot continuum emission of RU Lup at >1200 A is ap- 
proximately fitted by a blackbody with a temperature of 
6.8-7.8x 10 3 K and a luminosity of 2.7-6 L Q (Lamzin et al. 



1 We assume that the magnetospheric accretion columns do 
not significantly contribute to the radio flux due to its large 
optical thickness and small size, or obscuration by an optically 
thick jet/wind. The latter explanation is often used for the fact 
that WTTSs often show nonthermal radio emission from the 
magnetosphere, while it is usually absent in CTTSs and Class 
I YSOs (see, e.g., Andre et al. 1992; Martin 1996). 



1996), and a flux of ionizing photons of 10 37 -10 38 s _1 from 
this region could sustain only 10 _4 -10 _3 % of the recom- 
bination rate corresponding to the observed Ha luminos- 
ity. Thus, one of the remaining two mechanisms, X-ray or 
MHD-wave heating, would be more likely for the heating 
mechanism in this region. The same mechanisms are also 
proposed for magnetospheric accretion columns, which are 
responsible for Balmer and other permitted lines in many 
PMS stars (Lamzin et al. 1996; Calvet & Gullbring 1998). 

From the spatial scales, flow velocities, and Ha lu- 
minosities, we can estimate the typical hydrogen number 
densities and opening angles in these AU-scale outflows. 
The details of the calculation is described in Appendix 
B. Throughout we assume the following: (1) an elec- 
tron temperature of 10 4 K, (2) viewing angle of 45°, (3) 
M out = 0.1 x M acc , where M out and M acc are the mass 
loss rate and mass accretion rate, respectively, and (4) 
the flow has axisymmctric geometry and zero width at 
its origin. The first assumption is based on models which 
include X-ray heating (Shang et al. 2002) or the tempera- 
ture of magnetospheric accretion columns (e.g., Muzerolle 
et al. 2001), in which the presence of MHD-wave heat- 
ing is suggested by Calvet & Gullbring (1998). The sec- 
ond assumption is likely to be a good approximation for 
order-estimates, since viewing angles close to pole-on and 
edge-on would not provide the observed displacements and 
velocities, respectively. The third assumption is based on 
Calvet (1997), who revised the estimate from forbidden 
line and accretion luminosities by Hartigan et al. (1995). 
For mass accretion rates in RU Lup and CS Cha, we adopt 
7.7xl0 7 and 1.6xl0 7 M Q yr _1 , respectively, obtained by 
John-Krull, Valenti, & Linsky (2000). 

Fig. [To] shows the derived hydrogen number densities 
and opening angles in RU Lup and CS Cha, as a function 
of the assumed ionization fraction. The upper figure shows 
that the hydrogen number density is <; 10 cm -3 for any 
case, and larger for lower ionization fractions. Previous 
imaging observations have shown that opening angles of 
jets are likely to be <;20° within ~100 AU of the stars 
(see Eisoffel et al. 2000, and references therein; see also 
Bacciotti et al. 2000), and the lower figure shows that this 
fact gives lower limits to the ionization fraction of 0.2-0.5 
for the assumed volume filling factors of 0.1-1, providing 
the upper limits of n#^10 9 cm~ 3 . 

5.5. Disc gap/hole in CS Cha and RU Lup — 
implications for planet formation 

Recent radial velocity surveys have revealed a number 
of extrasolar giant planets (see, e.g., Marcy, Cochran, & 
Mayor 2000). These observations show a population min- 
imum between stellar-substellar companions and extraso- 
lar giant planets ("the brown dwarf desert" — Marcy, 
Cochran, & Mayor 2000; see also Halbwachs et al. 2000), 
suggesting the presence of another formation mechanism 
for the planetary-mass objects. The proposed mechanisms 
include: (1) core accretion via accumulation of dust grains, 
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ionization fraction 

Fig. 10. Mean hydrogen density and flow opening angle 
as a function of the ionization fraction. The estimated 
opening factor depends on the volume filling factor in the 
flow, which are assumed to be 0.1 and 1. The hached area 
indicates unlikely opening angles (see the text). 



which is generally accepted and requires 10 6 -10 7 yrs for 
planet formation (e.g., Pollack ct al. 1996), and (2) disc 
instability, which could form a gas giant planet only in a 
few hundred years (Boss 2000). 

Model calculations have shown that a young planet 
could tidally interact with a circumstellar disc, opening up 
a gap or hole at its orbital scale (see Lubow & Artymowicz 
2000; Takcuchi, Miyama, & Lin 1996). Even the most re- 
cent high-resolution techniques cannot resolve such AU- 
scale structures in the nearest star- forming regions, while 
infrared spectral energy distributions (SEDs) in several 
young stars exhibit dips at the mid-infrared, suggesting 
the presence of such gaps/holes (e.g., Strom et al. 1989; 
Marsh & Mahoney 1992,1993). Recent studies by Boss & 
Yorkc (1993,1996) show that the SEDs cannot provide 
convincing proof for the presence of a disc gap or hole 
as these dips could also be caused by a combination of 
vertical temperature structure and changes in dust opac- 
ity. 

On the other hand, our results for RU Lup (Takami et 
al. 2001) and CS Cha (this work) are providing evidence 
for the presence of disc gaps/holes at AU scales. These ob- 
jects exhibit the presence of Ha bipolar jets at AU scales, 
and the detection of the redshifted flows is explained by a 
disc gap/hole allowing this flow to be seen (Takami et al. 
2001; §3.4). Our position spectra would also be capable 
of discovering a stellar companion. No evidence for com- 



panion is seen, suggesting that the gap and hole in these 
objects could be induced by unseen companions. 

The SEDs suggest that the outer radii of the gap 
and hole correspond to a temperature of ^200 K, coin- 
cident with the ice condensation temperature (160-200 
K — see Marsh & Mahoney 1992; Boss 1995). Marsh 
& Mahoney (1992,1993) also obtained similar tempera- 
tures from mid-infrared dips in GK Tau and HK Tau 
SEDs for which stellar companions have not been de- 
tected. In contrast, the gaps/holes induced by a stellar 
companion show a wide range of temperatures: Guilloteau, 
Dutrey, & Simon (1999) provide an inner temperature of 
the GG Tau ring of ~ 30 K, while results by Jensen & 
Mathicu (1997) suggest temperatures up to ^600 K in 
the gaps/holes towards spectroscopic binaries. Marsh & 
Mahoney (1992,1993) suggest that such results could be 
explained if a gas-giant planet formed at the ice conden- 
sation radius and clearing up the innerside of the circum- 
stellar disc by tidal interaction. Indeed, theoretical work 
shows that the formation time scales of the gas-giant plan- 
ets are likely to be shorter at the ice condensation radius 
or outer region than the inner region: the presence of water 
ice could enhance the surface density of the solid material 
by a factor of 3-4 (e.g., Hayashi 1981), providing a shorter 
timescale for the planetesimal accretion (Lissauer 1987). 
Furthermore, water vapor produced in the inner region 
could condensate at the ice condensation radius, enhanc- 
ing the surface density of ice by a factor of ^70 (Stevenson 
& Lunine 1988), thereby allowing the ice condensation ra- 
dius to be a preferential location for formation of giant 
planets. 

RU Lup, CS Cha, GK Tau, and HK Tau have stel- 
lar ages of l-3xl0 6 yr (see Palla & Stahler 1999 for the 
evolutionary tracks, and Lamzin et al. 1996, Gauvin & 
Strom 1992, and Marsh & Mahoney 1993 for their effective 
temperatures and luminosities), implying that the young 
planets could be formed within this timescale. This re- 
sult contrasts with the detailed simulations by Pollack et 
al. (1996), who suggest that ^5xl0 6 yr is required for 
the formation of a Jupiter in the solar nebula. This age 
problem could be overcome if the initial surface density 
is higher than their standard model (10 g cm~ 2 for solid 
material at the scale of the Jovian orbit). These authors 
show that the timescale for the formation of giant plan- 
ets appears to be extremely sensitive to the initial surface 
density of the protoplanetary disc: e.g., increasing the den- 
sity by a factor of 2 could reduce the timescale by a factor 
of 30 for the formation of Jupiter. Such higher initial sur- 
face densities are likely to provide higher planetary mass 
than Jupiter (Lissauer 1995; Pollack et al. 1996), consis- 
tent with the discovery of a number of extrasolar plan- 
ets with such masses (see Marcy & Butler 2000; Jorissen, 
Mayor, & Udry 2001). We thus conclude that the disc gaps 
and holes observed in RU Lup and CS Cha are consistent 
with the scenario that these are induced by a gas-giant 
planet formed at the ice condensation radius. 
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6. Conclusions 

We present spectro-astrometric observations for 28 south- 
ern PMS stars. The structures detected in the position 
spectra include: (1) almost all the known binary compan- 
ions in our sample (Sz 68, Sz 41, HO Lup, VW Cha, S 
CrA, and AS 205), (2) companion candidates which have 
not been detected by infrared speckle techniques (T CrA, 
MWC 300), (3) monopolar and bipolar jets (AS 353A, CS 
Cha), (4) a combination of a jet and stellar companion (R 
CrA), and (5) a combination of jets and a bow shock (VV 
CrA). 

We succeeded to detect almost all the known binaries 
in our sample (6 out of 7) with separations of 0.1-1.5 arc- 
sees. The position of Ha emission in these objects are dis- 
placed primarily in the direction of the faint component, 
and their angular scales are consistent with those obtained 
by Bailey (1998b): 0.07-0.7 times as large as the binary 
separation. The displacement in the other features is much 
shorter than those in the Ha emission (0.02-0.05 times as 
large as binary separation) due to larger contribution of 
the continuum flux. All the observed displacements lie in 
the direction of the binary system, and these are explained 
by different spectra between the two binary components. 
Among the "spectro-astrometric" binary candidates, we 
confirm the time variation of the position angle in T CrA, 
consistent with binary motion. 

From results in known binaries to date, we conclude 
that this technique should be able to detect binaries with 
separations down to ~10 milliarcsec. This detection limit 
is comparable to that of lunar occultations, although our 
technique can easily achieve better limits by observing a 
larger number of photons (§2), and is applicable for objects 
with any coordinate and binary position angles. The Ha 
profiles obtained for the components of the known binaries 
are strikingly similar to each other — similar lincwidths 
and absorption features at similar velocities, results which 
support the scenario of core fragmentation for the mech- 
anism of binary formation. 

The Ha outflows observed by Bailey (1998b), Takami 
et al. (2001), and this work exhibit velocities coincident 
with the high-velocity components of the forbidden lines, 
suggesting that these are associated with jets close to the 
source. Among these objects, Z CMa and AS 353A show 
spatial scales more than 50 AU, and these are likely to be 
excited by shocks in the internal working surfaces. On the 
other hand, bipolar outflows in RU Lup and CS Cha ex- 
hibit much shorter spatial scales (1-5 AU), and the spatial 
scale and line profile in these objects suggest the presence 
of a heating mechanism other than shocks: e.g., X-ray or 
MHD-wave heating. From the observed spatial scales, flow 
velocities, and Ha luminosities of these flows, we estimate 
typical hydrogen number densities <; 10 7 cm~ 3 . 

Our results in RU Lup (Takami et al. 2001) and CS 
Cha (this work) suggest the presence of disc gaps/holes at 
AU scales, consistent with the presence of a mid-infrared 
dip in their SEDs. The SEDs suggest that the outer radii 
of these gap and hole have temperatures of ^200 K, co- 



incident with the ice condensation temperature (160-200 
K). The constant temperature in these objects and a few 
more disc gap/hole candidates contrasts with the wide 
range for those in PMS binaries (30-600 K). The presence 
of such disc gaps/holes could be explained if a gas-giant 
planet formed at the ice condensation radius, clearing up 
the innerside of the circumstcllar disc by tidal interaction. 
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Table A.l. Estimated radio luminosities 
object ifl El 



(mag) (A) 



RU Lup 
CS Cha 



9.99 
10.84 



29.7 
28.3 



A v 
(mag) 

1.28 
0.85 



(Lo) 



1.4xl0" 2 
3.6xl0 -3 



(mjy kpc 2 



0.42 
0.11 



a Herbig & Bell (1988) 

b sum for the blueshifted and redshifted emission in 
which positional displacement is observed 

c adopting the distances obtained by Wichmann et al. 
(1998) 

d Hughes et al. (1994) 

e Gauvin & Strom (1992) 



Appendix A: Expected luminosities for thermal 
radio emission 

The fluxes of Ha and free-free emission per unit volume 
are described as follows (see Osterbrock 1989; Reynolds 
1986): 

^]Ha = 3.5 x 10" 25 n p n e erg cm" 3 s" 1 (A.l) 
Anj v = 3.3 x 10" 38 n p n e i/" 01 erg cm" 3 s" 1 Hz" 1 ^) 

where n p and n e are the proton and electron densities, re- 
spectively; v is the observed frequency for free-free emis- 
sion. Throughout we assume an electron temperature 
of 10 4 K, plausible for either shock heating (see, e.g., 
Hartigan, Morse, & Raymond 1994b), ambipolar diffu- 
sion (Saher 1993; Garcia et al. 2001), or UV photoion- 
ization (e.g., Osterbrock 1989) as a heating mechanism. 
We also assume that the Ha emission is produced via 
recombination, and the contribution from collisional ex- 
citation is negligible. Indeed, the transition level of Ha 
is too high to be excited thermally at the assumed tem- 
perature (Osterbrock 1989). Collisional excitation of Ha 
could occur in fast shocks, although the total Ha flux in 
these shocks may still be dominated by the "recombina- 
tion plateau", which has a temperature of ~10 4 K and 
a much larger volume (e.g., by a factor of 10 3 ) than the 
region in which the line is collisionally excited (see Figs 
2-3 of Bacciotti & Eisloffel 1999). 

Equations (Al) and (A2) show that free- free emission 
is proportional to that of Ha emission, 



J3.6c 



= 9.4 x 10" 15 



JHo 



Hz" 



(A.3) 



Table A.l summarizes the observed parameters, and the 
estimated luminosities for free-free emission for the opti- 
cally thin case. Note that the derived free- free luminosities 
are upper limits, since the radio thermal emission in YSOs 
are often optically thick (see Beltran et al. 2001, and ref- 
erences therein). 



Appendix B: Hydrogen density and opening angle 
of the Ha bipolar outflow 

The total mass of the Ha emission line region and the Ha 
luminosity are described as follows: 



m p n H dV = M out l out /v out 



(B.l) 
(B.2) 



LHa = 3 ' 5Xl(r25 / npUedV CTSS " 
where m p and nu are the proton mass and the hydro- 
gen number density, respectively; M out is the mass loss 
rate; l ou t and v ou t are the spatial scale and the velocity of 
the observed outflow; n p and n e are the proton and elec- 
tron number densities, respectively, in units of cm" 3 ; V 
is the volume of the Ha emission region. Equation (B2) 
is derived adopting the recombination rate in Osterbrock 
(1989) assuming an electron temperature of 10 4 K (see 
Appendix A). We also assume that collisional excita- 
tion of Ha is negligible, since its transition level is too 
high to be thermally excited at the assumed temperature 
(Osterbrock 1989). Using equations (Bl) and (B2), the 
mean electron density is derived as follows: 



(n e ) 



J n p n e dV _ L Ha [erg s x ] m p v out 



Jn p dV 3.5 x 10" 25 M out l out ( Xi ) 

I \-i ( Llia \ ( v ° ut \ 
{Xi) I 1 L m J Vl00kms-V 



1.9 x 10 



10 



-l 



10" 8 M yr" 1 



1 AU 



cm" 3 (B.3) 



where the mean ionization fraction is defined as = 
(J n p dV)/(J nu dV). From the equation (B3), the mean 
hydrogen number density is derived as follows: 

(n H ) = (n e )/(xi) 



1.9 x 10 10 ( Xi y 



1 L 



Vout 



10" 8 Mq yr" 



100 km s- 1 

^out 



1 AU 



-3 



(B.4) 



Next, we derive the typical opening angle of the flow as- 
suming axisymmetric geometry and that the flow has zero 
width at its origin. The total mass in the flow is described 
as follows: 



M out ~ / Kr 2 out m p n H f H dl 



^(r 2 out )l out m p (n H ){f H ) 



(B.5) 



where r out is the radius (or half-width) of the flow at a 



certain distance from the star; 



t ) is the mean square 



of the radius; is the volume filling factor. On the other 
hand, the typical flow width (d out ) and opening angle (9) 
can be described as follows: 



{d ut) 

(0) e 2 Ian 



_i {d ut) 
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(B.6) 
(B.7) 
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Substituting equations (B1)(B5)(B6) to (B7), (9) is de- 
rived as follows: 



(6) ~ 2 tan" 1 ,/- M ° ut ,„ — , (B.8) 



